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ABSTRACT 

The available information on isotopic abundances in the atmospheres of low- mass Red Giant Branch 
(RGB) and Asymptotic Giant Branch (AGB) stars requires that episodes of extensive mixing occur 
below the convective envelope, reaching down to layers close to the hydrogen burning shell (Cool 
Bottom Processing). Recently iBusso et al.1 (|2007[ ) suggested that dynamo-produced buoyant mag- 
netic flux tubes could provide the necessary physical mechanisms and also supply sufficient transport 
rates. Here, we present an a — 51 dynamo in the envelope of an RGB/ AGB star in which shear and 
rotation drain via turbulent dissipation and Poynting flux. In this context, if the dynamo is to sustain 
throughout either phase, convection must resupply shear. Under this condition, volume-averaged, 
peak toroidal field strengths of (B^) ~ 3 x 10 3 G (RGB) and (B^) ~ 5 x 10 3 G (AGB) are pos- 
sible at the base of the convection zone. If the magnetic fields are concentrated in flux tubes, the 
corresponding field strengths are comparable to those required by Cool Bottom Processing. 
Subject headings: MHD - stars: AGB and post-AGB - stars: interiors - stars: abundances -stars: 
rotation - giant stars 



1. INTRODUCTION 

In low-mass Red Giant Branch (RGB) and Asymp- 
totic Giant Branch (AGB) stars, matter must circu- 
late in episodes of extended mixing from the convec- 
tive zone to the radiative region where nuclear processing 
can occur (the so-called Cool Bottom Processing; CBP; 
Wasserburg et alj|1995t ICharbonnel and Do Nascimentol 



19981: INollett et all 120031 : iHerwid 12001) . While the 
mechanism driving mixing remains under investigation, 
CBP predicts the chemical and isotopic evolution for 
a given bulk transport rate by utilizing the path in- 
tegral of the nuclear reactions over t he mass and in- 
dependent of the mixing mec hanism (jBoothrovd et all 
119941 : IWasserburg et all I1995T ) . The necessary trans- 
port rates and processing temperatures are constrained 
by RGB/ AGB isotopic measurements and presolar me- 
teoritic dust grain abundance ra t ios (IChoi et al. 119981 : 
lAmari et all l200l INollett et al.1 120031: iNittlerl 120051) . 
Weak circulation (M ~ 10~ 6 - 10~ 8 Af /yr) in which 
material is transported from the convective envelope to 
just above the H-burning shell, processed and then re- 
turned has had success in matching: (i.) 12 C/ 13 C, C and 
N abundances in RGB stars, (ii.) 12 C/ 13 C and N/O in 
AGB stars and (iii.) 26 A1/ 27 A1, 18 0/ 16 0, 17 0/ 16 and 
12 C/ 13 C in circumstellar dust grains f ound in meteorites 
that must come from AGB sources (IWasserburg et al.l 
[19951 INollett et al.1l2003l : lHerwidl2005h . 

Several physical origins have been proposed for 
the assumed mixing including: shear and thermo- 
haline instabilities, meridiona l circulation and inter- 
nal gravity waves (jZahnl 119921 : iDenissenkov and Weiss! 
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1 19961: [De nissenk ov and Tout! 120031: IChaname et all 
120051 lEggleton et al.ll2006l ; ICharbonnel and Zahnll2007h . 

Purely rotation-induced mixing appears insufficient as 
isotopic changes are minimal during the RGB and lead 
to a quenching of s-process nucleosynthesi s during the 
AGB (iSiess et al.ll200l iPalacios et al] 12006ft . 

Recently, it was proposed that buoyant magnetic flux 
tubes can induce the necessar y transport rates r equired 
for Cool Bottom Processing ([Busso et alj|2007ft . Con- 
straints on the field strengths for the required mixing 
rates were calculated while the fields were assumed to 
originate from a dynamo operating in the RGB and AGB 
interiors. In particular, the proposed mechanism would 
provide a means of rapid transport of material to the con- 
vective mantle and may provide a solution to the problem 
of Li rich stars. 

Magnetically mediated outflows resulting from dy- 
namo amplification have been proposed as the origin 
of bipolarity in post-Asympotitc Giant Branch stars 
(post-AGB) and plane t ary nebula fPN) dPascolil 1 19971 : 
Blackm an ct al] l200ll iNordhaus and Blackmanf 120061 : 
Nordhaus et al.|[2007f ). However, dynamo-induced shap- 
ing requires differential rotation. This can be supplied 
by binary companions or sustained throughout the AGB 
phase by convective re distribution of angular momentum 
(|Nordhaus et al.ll2007h . Convection may resupply differ- 
ential rotation (analogous to the A-effect in the Sun) and 
sustain the magnetic fields throughou t the RGB /AGB 
stage s ([Riidiger and Hollerbachl 120041 : INordhaus et al.1 
2007). In particular, it may be that a weaker dy- 
namo (supplying sufficient mixing) is sustained through- 
out the RGB/ AGB phases, but that a binary compan- 
ion is required to power bipola r jets in post-AGB stars 
(|Nordhaus and Blackmanll2008fh 

We investigate the origin of large-scale magnetic fields 
via a dynamical, a — 51 dynamo operating at the base of 
the convection zone in an initially 1.5 Mq main sequence 
star during the RGB and AGB phases. Our s tellar mod- 
els are identical to those used in IBusso et al.l (|2007l ) , and 
allow for direct comparison. In both phases, subadia- 
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batic zones similar to the solar tachocline exist. We al- 
low a fraction of the turbulent energy cascade to resupply 
shear, thus sustaining the dynamo. The field penetrates 
the shear zone and may transport material close to the 
H-burning shell. We compare our results to the field 
strengths and penetration depths needed for adequate 
magnetic mixing in low-mass RGB and AGB stars. 

2. MAGNETIC MODEL 

A 2-D schematic of our 1-D model is presented in Fig. 
[TJ The convection zone extends from the stellar sur- 
face to the interface between the convective and radia- 
tive zones. Convective twisting motions convert buoyant 
toroidal fields into poloidal fields through the a-effect. 
Below the convection layer, the differential rotation zone 
shears poloidal fields back into toroidal fields via the fi- 
effect. 
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Fig. 1. — Geometric representation of the stellar interior. The 
poloidal field is amplified in the convective zone, while toroidal 
field is amplified in the shear layer. The bottom of the convection 
zone is given by r = r c while the maximum depth required by Cool 
Bottom Processing is r = r c i, p . The top of the H-burning shell is 
located at r = r/,. 



We solve for the time evolution of dynamic quantities 
at the interface (r = r c ). The rotation profile across 
the shear layer varies from at the interface to f2 + 
Af2 at r = r c — L. Thus, Af2 is a measure of shear in 
the differential rotation zone. If Ail = 0, the system 
exhibits solid body rotation. The average poloidal field 
(Bp), and average toroidal field (B^), amplify from 1 G 
seed values until they are quenched through a drain of 
the avail able differential rotatio n energy. We refer the 
reader to iNordhaus et al.1 ((2007) for a derivation of the 
precise mean-field equations solved. 

2.1. Shear Zone Penetration Depth 

To capture aspects of the 2-D, mean-field geometry 
within the framework of our 1-D time-dependent model, 
we employ two turbulent diffusion coefficients: (3 P , cor- 
responding to diffusion of the poloidal field (which grows 
primarily in the convective region) and fi$, correspond- 
ing to diffusion of the toroidal field (which is amplified 
in the differential rotation zone; see Fig. 1). We also 
employ as the turbulent diffusion coefficient for the 
toroidal velocity. 

The convective region is highly turbulent and the dif- 
ferential rotation zone is weakly turbulent, therefore 



(3<f> <C P p . The value of /3$ determines how far the toroidal 
component of the field can diffuse into the shear zone in a 
cycle period. The further into the shear zone the toroidal 
field can penetrate, the greater the shear energy that can 
be extracted and utilized by the dynamo. 

The depth to which the toroidal field can diffuse into 
the shear layer in a cycle period, r, is defined as S ~ 
(A/>t) 5 . The cycle period does increase in the dynamical 
regime, however, it does not change appreciably from its 
initial value. We set 5 — L so that the field diffuses to the 
depth required by Cool Bottom Processing (r — r c bp). 
This fixes (3^, : rgb — 5.6 x 10 14 cm 2 s _1 and P^.agb — 
6.9 x 10 15 cm 2 s _1 as the required transport rates are 
~ 5 flux tubes per year during the RGB and ~ 8 flux 
tubes per year during the AGB (jBusso et al.ll2007l ). 

2.2. The Need For Convection Resupplying Shear 

As the magnetic field is amplified, differential rota- 
tion energy is drained to support its growth. The 
back-reaction of field amplification on the differential 
rotation results in rapid termination of the dynamo 
(< 100 yrs) from the initial available shear energy; see 
INordhaus et al.l ((20071 ) for details. In order to sustain 
the dynamo through an RGB/ AGB lifetime, a constant 
differential rotation profile must be established. This 
occurs in th e sun as convection re-seeds s hear through 
the A-effect ((Riidiger and Hollerbachll2004( ). Although it 
remains to be established if a similar effect occurs in 
evolved stars, by analogy to the solar case, we allow 
a fraction of the turbulent energy cascade to resupply 
shear. Additionally, we keep the rotation at the interface 
fixed. This is physically equivalent to storing the Poynt- 
ing flux in the interface region. If the field is trapped, the 
Poynting flux does not emerge from the layer and thus, 
does not spin down the envelope. An RGB/ AGB dynamo 
will remain stable when the following two conditions are 
met: (i.) convection resupplies s hear (ii.) the Poyntin g 
flux is stored inside the envelope (|Nordhaus et al.ll2007h . 

The differential rotation profile in the interior is un- 
known. An estimate for the maximum sustainable shear 
that convection can resupply for a sufficiently rotating 

1/2 

star is given by AVt max = (M c /Mah) i v /L) where 
M c is the mass of the convective layer, M^n is the mass 
of the shear zone and v is a typical convective veloc- 



ity. For the RGB phase, A£l r < 



2.5 x 10" 



whilc 



A£l max ~ 5.5 x 10 -5 s _1 for the AGB phase. If the initial 
differential rotation were greater than AQ max , the mag- 
netic field would extract the excess shear before relaxing 
to sustainable values. In this case, an initial period of 
enhanced magnetic activity would be expected followed 
by steady dynamo action. In §3, for suitable Af2, we de- 
termine the fraction of turbulent cascade energy per unit 
time necessary to maintain a steady-state dynamo i n our 
RGB/ AGB model stars (e.g. INordhaus etal][2007l ). 

2.3. Flux Tubes vs. Average Fields 

Our calculations are for mean-fields and correspond 

to volume averages. When the ratio of average ther- 

(p) 

mal to magnetic pressure p = rgW7§^ > 1 the mag- 
netic field may be concentrated in flux tubes with an 
overall volume filling fraction £, such that ~ 4 

l|Blackmanlll996[ ). When the magnetic pressure in flux 
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TABLE 1 

CONVECTIVE AND SUBCONVECTIVE PARAMETERS 



Parameter 


Symbol 


RGB 


AGB 


Base of conv. zone (cm) 


r c 


6.3 x 10 10 


5.4 x 10 10 


CBP radius (cm) 


^cbp 


3.4 x 10 9 


1.9 x 10 9 


H-burning radius (cm) 


Th 


2.0 x 10 9 


1.5 x 10 9 


Mass of conv. zone (g) 


M c 


2.2 x 10 33 


1.1 x 10 33 


Mass of C-layer (g) 




1.0 x 10 31 


1.3 x 10 30 


Rot. at interface (s _1 ) 


n(r c ) 


3.0 x 10" 6 


2.5 x 10~ 6 


Diff. Rotation (s^ 1 ) 


An 


1.01 x 10~ 5 


1.67 x 10- 5 


Model Results 


Toroidal Field (G) 




2602 


4778 


Poloidal Field (G) 


(B P ) 


342 


596 


Re-supply Rate 


f 


0.154 


0.075 


Vol. Filling Fraction 


c 


0.36 


0.23 


CBP Requirements" 


Flux Tube Field (G) 


(Bt,c) 


3450 


8600 


Flux Tube Field (G) 


( B t,cb v ) 


3.8 x 10 5 


5 x 10 6 



"From lBusso et alj IpOO^I 



tubes are in pressure balance with their exterior, the 
mean-field strength is then related to a flux tube field, B t 
by (B 2 ) ~ (Bt)- Sinc e the requirements from mag- 
netic mixing (jBusso et al.ll2007h constrain the magnitude 
of flux tube fields to have C < 1, the mean field can be 
lower than the flux tube field, relaxing the demands on 
the strength of the dynamo and the required shear. 

3. DYNAMO AMPLIFIED MAGNETIC FIELDS 

By requiring CBP transport rates, the necessary field 
strengths and buoyant velocities are calculated for a 
given stellar mod el. We employ the same models as 
iBusso et~a l. (2007), in which an initial 1.5 M Q main se- 
quence progenitor is evolved through the RGB and AGB. 
For each phase, we determine what is required to sustain 
the dynamo and the corresponding saturated large-scale 
field strengths from our a — f2 dynamo. A summary of 
the model parameters can be found in Table [T] 

3.1. The RGB Case 

For our model RGB star, the stellar mass is relatively 
unchanged from the main sequence [M = 1.49 M©). The 
base of the convection zone is located at r c = 6.3 x 10 10 
cm while the maximum depth to which matter must cir- 
culate is r c b p — 3.4 x 10 9 cm. This implies a shear 
zone thickness of L ~ 6 x 10 10 cm. The mass of the 
shear layer and convective envelope for the RGB star are 
M An = 1.0 x 10 31 g and M c = 2.2 x 10 33 respectively. 
The rotation rate at r = r c is obtained by conserving 
angular momentum on spherical shells from a main se- 
quence pro genitor whose surface rotation velocity is ~ 
30 km s" 1 (|Kawaler|[l987T) . 

Using CBP mass transfer rates, mixing depths and 
buoyant rise velocities, the necessary magnetic field 
strength at r c and r cbp can be constrained. If the field 
is purely toroidal and in flux tubes that rise individu- 
ally, then the average toroidal flux tube field required 
by CBP is (B t , cbp ) = (B t (r cbp )) = 3.8 x 10 5 G. At the 
base of the convection zone (assuming conservation of 
mass and magnetic flux), this corresponds to a field of 



B t .c) = (B t (r c )) ~ 3.5 x 10 3 G for a dipole geometry 
Busso et ai1l2007| ). 

To compare, we calculate the mean field at the base 
of the convection zone assuming that a fraction /, 
of the turbulent ener gy cascade rate resupplies shear 
(jNordhaus et aLlfeOOTh . The value of / is not known, 
nor is it established that convection resupplies shear in 
evolved stars. Although our choice of / is arbitrary, we 
purposely pick low values so that the requirements on 
the energy that convection must supply are reasonable. 

For the parameters in TableQ] the peak toroidal field at 
the base of the convection zone is (B^) = 2.6 x 10 3 G (see 
Fig. [2]). In the RGB star, CBP requires approximately 
5 buoyant flux tubes per year to reach the convective 
envelope. This is comparable to the dynamo cycle period 
of ~ 0.2 yrs and suggests the necessary transport rates 
are plausible. To sustain the dynamo throughout the 
RGB lifetime, ~ 15% of the turbulent cascade energy 
must reinforce shear. This is not a lower limit as smaller 
resupply rates will also sustain the dynamo, albeit with 
lower peak field strengths. 

To compare our mean-fields with the flux tube fields 
necessary for CBP, we determine the volume filling frac- 
tion, C, as detailed in §2.3. For the RGB model, we 
find that C = 0.36. However, it should be noted that 
the convective resupply rate is the dominant quantity 
as it sets the differential rotation and hence the peak 
field strengths. A smaller /, implies a lower (B^) which 
requires a smaller Q to match the required flux tube 
strengths. Further work on how convection resupplies 
shear in evolved stars is warranted. 
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Fig. 2. — Toroidal field as a function of time at the base of the 
convection zone for our model RGB star. The field amplifies until 
it saturates. 



3.2. The AGB Case 

The total mass of our AGB star has decreased from the 
initial 1.5 M© progenitor to M — 1.2 M©. The base of 
the convection zone is located at r c = 5.7 x 10 10 cm while 
the maximum depth required for CBP is r cbp — 2 x 10 9 
cm. The shear zone thickness has slightly decreased from 
the RGB phase such that L ~ 5.2 x 10 10 cm. The mass 
of the convective layer and shear zone have decreased to 
M c = 1.1 x 10 33 g and M An = 1.3 x 10 30 g respectively. 
The rotation rate at the interface is given by f2 (r c ) = 
2.5 x 10~ 6 s _1 while the shear is slightly larger in the 
AGB phase. 

A strong toroidal field is required at the CBP radius 
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Time (yrs) 

Fig. 3. — Toroidal field at the base of the convection zone for 
our model AGB star. 

{B t ,cb P ) = 5 x 10 6 G such that ~ 80 flux tubes per 
year rise to the convective envelope. At the base of the 
convection zone, the required flux tube field strength is 
(B tjC ) ~ 8.6 x 10 3 G (IBusso et alJl2007f ). 

From our dynamo calculations, the peak toroidal field 
at the base of the convection zone is (B^) = 4.8 x 10 3 G. 
Figure [3] shows the time evolution of the mean toroidal 
field in our model AGB star. The cycle period is ~ 0.01 
yrs which is comparable to the required buoyant flux tube 
frequency. 

To sustain the dynamo, ~ 7% of the turbulent energy 
cascade must rcsupply differential rotation. In general, 
more convective energy is available in the AGB phase, 
and thus, the resupply rate necessary to support differ- 
ential rotation is lower compared to the RGB star. 

Our mean fields require £ = 0.23 to match CBP flux 
tubes at the base of the convective zone. However, we 
note again that a lower /, will produce weaker peak field 
strengths and, in turn, a lower volume filling fraction. 

4. CONCLUSIONS 

IBusso et al.l ()2007j ) recently suggested that buoyant 
toroidal fields, may supply the necessary chemical trans- 
port rates to circulate material from near the hydrogen 
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burning shell to the convective envelope required by Cool 
Bottom Processing in low-mass RGB and AGB stars. 

In this paper, we have shown how an a — dynamo op- 
erating in the interior of an RGB /AGB star could supply 
the fields necessary to accomplish the needed transport. 
The fields diffuse in a magnetic cycle period to the depths 
required by Cool Bottom Processing. If the dynamo is to 
sustain throughout either phase, convection must resup- 
ply differential rotation. This can be accomplished if a 
fraction of the turbulent energy cascade reinforces shear 
(< 15% for the RGB star; < 7% for the AGB star). 

For the model RGB star, mean toroidal field strengths 
of (B$) = 2.6 x 10 3 G are sustainable at the base of the 
convection zone. If the magnetic field is concentrated 
in flux tubes, the volume averaged fields, (-B 2 ), are re- 
lated to the flux tube fields, {B 2 ), by Q < 1, such that 
(B 2 ) ~ (B 2 ). For the model RGB star, ( = 0.36 
will satisfy GBP requirements. For the model AGB star, 
stronger fields are possible with sustained toroidal field 
strengths of (B$) = 4.8 x 10 3 G possible at the base of 
the convective zone. This requires £ = 0.23 to match 
CBP flux tube requirements. 

We have found that dynamo-produced fields provide 
a plausible mechanism for the physical origin of extra 
mixing in low-mass RGB and AGB stars. However fu- 
ture work, in particular, realistic rotation profiles, the 
complex relationship between convection and shear, in- 
termittent magnetic mixing, the dynamics of mass trans- 
port via flux tubes all need to be studied further. Equally 
important is to independently constrain the value of the 
diffusion coefficients in the shear and convection zones to 
determine if they are consistent with what our model im- 
plies. Finally, > 2-D studies of this problem are needed. 
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